Pair Decay Width of the Hoyle State and its Role for Stellar Carbon Production 
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The pair decay width of the first excited + state in 12 C (the Hoyle state) is deduced from a novel 
analysis of the world data on inelastic electron scattering covering a wide momentum transfer range, 
thereby resolving previous discrepancies. The extracted value T w = (62.3±2.0) peV is independently 
confirmed by new data at low momentum transfers measured at the S-DALINAC and reduces the 
uncertainty of the literature values by more than a factor of three. A precise knowledge of rV 
is mandatory for quantitative studies of some key issues in the modeling of supernovae and of 
asymptotic giant branch stars, the most likely site of the slow-neutron nucleosynthesis process. 
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X 

=5 



(N 
> 

oo 
m 

o 
o 



The production of the most abundant stable carbon 
isotope 12 C in the center of stars proceeds through the 
fusion of three a particles (the triple-a process) . Its re- 
action rate is of critical significance for a variety of key 
problems of nuclear astrophysics jl| like the elemental 
abundance in the universe Q , the size of the iron core in 
massive stars determining the properties of supernova ex- 
plosions 0, the dynamics of asymptotic giant branch 
(AGB) stars [5[, the site of the main slow- neutron cap- 
ture process (s- process) of heavy elements @ , or the weak 
s-proccss in massive stars 0. 

Under most astrophysical conditions the reaction ex- 
clusively proceeds through a scattering resonance of three 
a-particlcs that represents an excited state in 12 C with 
quantum numbers J*" = + at an excitation energy 
E x = 7.654 McV (the so-called Hoyle state 3), which 
then decays to the stable ground state. The reaction 
rate of the resonant triple-a process is proportional to 
the radiative decay width T ra d of the Hoyle state. Al- 
though this state is experimentally studied for more than 
50 years, at present T ra d is known with an uncertainty of 
about ±12% only. However, an accuracy of about ±5% 
is requested for applications in astrophysics 0, [B[ . 

Experimentally the radiative width cannot be accessed 
directly but is determined as a product of quantities 
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measured in different experiments. Here, T — T a + T 7 + 
is the total decay width taking into account a, 7 and 
e ± decays of the Hoyle state. Whereas the first quantity 
on the right-hand side of Eq. has been precisely de- 
termined with an uncertainty of ±2.7% (see Ref. |9| and 
refs. therein), the branching ratio T^/T is known with 
an error of ±9.2% only (see Ref. [l(| and refs. therein). 
However, a new measurement of T v /T with an exp ected 
precision of about ±5% is currently underway 



The pair decay width in Eq. ([T} can be determined 
by inelastic electron scattering through the relation [l2[ 
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e c 2 ) 2 , a de- 
notes the fine structure constant and (r 2 ) tr the ma- 
trix element of the monopole transition. The energy- 
dependent correction term B(E X ) is given in Tab. 1 of 
Ref. [12| (0.898 for the present case). The influence of the 
nuclear Coulomb field is taken into account by the factor 
C(Z,E X ) and amounts to 1.013 for the Hoyle state. The 
reduced £^0 transition probability is the largest known 
for excitation of a single state [13| and exhausts about 
7.5% of the energy-weighted sum rule. 

Two values were extracted from measurements at 
low momentum transfers q with a so-called 'model- 
independent' analysis [l4|; 1 15| in the distorted wave 
Born approximation (DWBA) explained in detail below. 
They agree well with each other and an averaged value 
60.5 ± 3.9 ^cV is quoted in Ref. [3]. Alternatively, a 
Fourier-Bessel analysis 17j of the transition form factor 
including data over a wide momentum transfer range has 
recently been done by Crannell et al. [3| ■ The extracted 
value of r w = 52.0 ± 1.4 ^eV has an error of 2.7% only, 
but the \ow-q result deviates by about 6tr. 

To resolve this discrepancy we have performed a new 
high-resolution electron scattering experiment at low mo- 
mentum transfer and reinvestigated both approaches. A 
novel model-independent ansatz was developed to ana- 
lyze the global form factor. This provides not only a 
precise determination of the pair width T w but also an 
important test of current models for the structure of the 
Hoyle state (lij . which represents a prime example of 
an a-cluster state predicted to possess the features of a 
low-density gas of a particles resembling a Bose-Einstein 
condensate |20j|. As shown in Ref. Il9f, state-of-the-art 
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models describe the (e,e') form factor quite well over a 
broad momentum transfer range but predict a monopolc 
matrix element (directly related to T^) about 20% too 
large (cf. Table I in |l9|). The present very precise re- 
sults confirm this discrepancy and provide new insight 
into its origin. 

The experiment was carried out at the high energy- 
resolution magnetic spectrometer [2l| of the Darmstadt 
superconducting electron linear accelerator S-DALINAC. 
Excitation energy spectra were taken at initial electron 
energies between 29 McV and 78 MeV and scattering an- 
gles from 69° to 141° with beam currents of about 1 p,A. 
The momentum transfer range for the excitation of the 
Hoyle state varied thus between 0.21 fm _1 and 0.66 fm _1 . 
A self-supporting carbon target (98.9% 12 C) with an 
areal density of 6.4 mg/cm 2 was used. In dispersion- 
matching mode an energy resolution AE w 28 keV (full 
width at half maximum, FWHM) was achieved. Exam- 
ples of spectra at a beam energy of 73 MeV are presented 
in Fig. [1] The peaks correspond to the elastic line, the 
2f state and the Hoyle state in 12 C. 
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FIG. 1: Spectra of the 12 C(e,e') reaction measured at a beam 
energy Eo = 73 MeV and scattering angles 9 = 93° (top) and 
6 = 141° (bottom). 

As the form factor is the Fourier transform of the 
density, it corresponds to a cross section calculated in 
plane wave Born approximation (PWBA). The measured 
cross section, however, corresponds to the distorted wave 
Born approximation (DWBA). Therefore, we deduce the 
square of the experimental form factor - the differential 
cross section divided by the Mott cross section - at the 
measured momentum transfer qi and electron beam en- 
ergy E oi by 
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where Z is the number of protons and B(C0,qi, E^i) 
denotes the reduced transition probability, while 
B PWBA ( ?i ) and B DWBA ( ft ,So i ) are reduced transition 
probabilities calculated within plane wave and distorted 
wave Born approximation, respectively. This so-called 



PWBA transformation allows to relate cross sections 
measured at different qi and Eoi directly to the form fac- 
tor or corresponding transition density. The data of the 
different experiments should then collapse into a single 
line for |.F t r xp (<7i)| . For the present measurements typical 
transformation factors are around 0.85. 

The factors B PWBA { qi )/ B DWBA ( qil E Ql ) in EajEJ are 
computed with the code of Heisenberg and Blok [22j in an 
iterative procedure starting with a the transition density 
Ptr{r) taken from the models discussed in Ref. flij. e.g . 
the Fermionic Molecular Dynamics (FMD) model |23| . 
The resulting form factor is transformed into an improved 
transition density that enters in a second step into the 
PWBA transformation in Eq. ([3]) etc. Convergence is 
reached after three steps. 

The transition form factor F tT (q) of a monopole tran- 
sition is related to the transition density ptr{r) by 
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where ptr(r) = (Of |/S(r)|0^~) is the matrix element of the 
charge density operator p{f) between the ground state 
and the Hoyle state and the operator Z = J p(r) d 3 r 
counts the number of protons. 

Expansion of the Bessel function jo(qr) for low q, i.e. 
<7 2 ( r2 ) tl . "C 1, shows that the form factor in Eq. Q 
is governed by the lowest moments (r 2 ) ti and (t 4 ) . 
Hence precise data at low q are important for an accu- 
rate determination of the pair width which is propor- 
tional to ( r 2 ) tr , cf. Eq. ([2]). But as we shall show later 
the power expansion around q = is very sensitive to 
statistical and systematic errors. Therefore we combine 
our new results with the previous world data set [l8| and 
perform a model-independent analysis of the 0^ — > 0^ 
transition form factor by a global fit to all data. 

In the present analysis we use the model-independent 
ansatz 
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which respects the condition F tr (q = 0) = and the 
fact that jo(qr) contains only even powers of q. The 
fit parameters are b and c n , where b plays the role of 
a length scale. The matrix element (r 2 \ = |6ci6 2 | is 
simply given by c\ and b. 

The transition density ptr(f) corresponding to the 
ansatz ((5|) is given by 
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where the relation between the coefficients c n and d n can 
be calculated by the inverse transformation of Eq. (jj]). 



For the Hoyle state measurements of |F t ° r xp (gi)| 2 exist up 
to (jWx = 3.1 fm -1 so that structures in pt T {r) can be 
resolved on a scale of ir/q max ~ 1 fm. 

Due to the Coulomb distortion of the in- and outgo- 
ing scattering states the cross section remains finite at 
momentum transfers where the form factor goes through 
zero and hence the uncertainty in deducing |F t cxp (<7)| ac- 
cording to the PWBA transformation Eq. §3§ is large 
around q = 2.1 fm -1 , see Fig. [5] Therefore only data 
F t e r xp (g)| 2 > 10~ 6 are considered. Variation of the max- 
imum power of the polynomial in Eq. ([5]) shows that 



10" 



10" 



10" 



- 10" 



10" 



10" 





1 1 1 1 

12 12 




C(e,e') C 




E x = 7.654 MeV 















q [fm" ] 

FIG. 2: (color online). Squared form factor for the transition 
from the ground state in 12 C to the Hoyle state extracted by 
the model-independent analysis (solid line) described in the 
text in comparison to the experimental data [TJ|, [7j| trans- 
formed according to Eq. ©. 

Figure [5] shows the resulting global fit. It describes 
the data well over the whole q range of measured mo- 
mentum transfers including the minimum region, where 
data were excluded from the fit. The following transi- 
tion radii are extracted: ( r 2 > fa = (5.47 ± 0.09) fm 2 and 
( r 4 V = (115 ± 8) fm 4 . The uncertainties are estimated 
by varying the data base because they are dominated 
by systematic errors in the data sets. The accuracy is 
essentially limited by the experimental uncertainties at 
q > 1.7 fm~ where the cross sections are small. This 
might look surprising because the transition radius is 
given by the limit q — > 0. However, we have reached 
an accuracy where seemingly small contributions become 
relevant and furthermore the form factor has to be a 
smooth function so that information from higher q values 
may influence the low-g behavior to a certain extent. 

In Fig. [3]we display pti(r)r 4 to show that ( r 2 V , which 
is 47T times the integral over this function, is very sensitive 
to contributions beyond r = 4 fm where the ground state 
density is very small but the density of the Hoyle state 
is not, see Ref. The transitions densities obtained 

in Fermionic Molecular Dynamics and BEC (gas of a- 
particles) are somewhat too large beyond 4 fm indicating 




FIG. 3: (color online). Transition charge density correspond- 
ing to the form factor displayed in Fig. [2]multiplied by r (full 
line) in comparison to the theoretical predictions [7J|] of FMD 
(dotted line) and an a particle gas (dashed-dotted line). 



that the calculated charge density of the Hoyle state is 
overestimated at large distance. Although the deviation 
docs not seem to be big the transition radius ( r 2 V ° c is 
for both models about 20% too large because negative 
and positive contributions cancel up to r w 3.5 fm and 
only the region beyond 3.5 fm matters. 

Crannell et al. 18[ used the same data, except for the 



seven new data points, and performed a Fourier Besscl 
(FB) analysis [17J. Figure |3] illustrates that beyond their 
adopted cut-off radius R c » 8.5 fm the transition density 
vanishes. They concluded however a significantly smaller 
r^. To understand this discrepancy we take our fit as 
an approximation of the experimental form factor and 
read off the FB expansion coefficients at q v = vn j R c for 
v = 1 up to v max = 8. The resulting transition density 
has (r 2 \ = 4.99 fm 2 or T™ = 51 ueV in accordance 
with the result of |18j . In order to reproduce our own 
result within a FB analysis we have to go up to at least 
^max = 10 which corresponds to q ma x = 3.7 fm -1 where 
no data exist. The same holds true if we repeat the anal- 
ysis without the S-DALINAC data. In general one should 
keep in mind that due to the oscillatory behavior of the 
Bessel functions a FB expansion can only reproduce the 
tail region well, if v max and q ma x are sufficiently large. 

The Taylor expansion of the form factor around q = 
offers itself a model-independent way to deduce the 
transition radius ^ r 2 V which is the leading coefficient 
in the expansion. To study this often employed approach 
we deduce from Eq. (T?]) 
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which is used to fit the low-g data measured at the S- 
DALINAC. The quantities <r 2 )|° W9 and (r 4 )|° w? are 
taken as free parameters. Higher powers up to q 8 are 
included in the fit such that the last term in Eq. 
contributes less than 1% of the first term. Because of 
the limited number of data points we approximated the 
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higher transition matrix elements by (( r 2 )J.° W<? )" ■ x n 
with x„ = ( r 2n V I ( r 2 )" . These ratios were taken from 
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FIG. 4: (color online). Experimental form factors at low 
q corrected by Eq. ((3j compared to different model ap- 
proaches. Extrapolation to zero momentum transfer provides 
the monopole matrix element. Square symbols and dashed 
line: new data and polynomial fit with Eq. (0. Circles and 
full line: old data and global fit with Eq. (J5j . Dashed-dotted 
line: FMD model, dotted line: a-cluster model fl{J. 

The dashed curve in Fig. [4] represents the fit to the 
data resulting in (r 2 ) 1 ™ 9 = (5.29 ± 0.14) fm 2 and 

( )'° W<? = (96 ± 10) fm 4 . As can be seen from Fig. 0] 
an extraction of the transition radius by extrapolation 
to q = considering only \ow-q data is very sensitive to 
the experimental errors. The global fit to all data (full 
line) reproduces the new low-g data except for a single 
data point at q 2 = 0.437 fm~ 2 . The slope of the poly- 
nomial fit (dashed line) is apparently not steep enough 
to match the data at higher momentum transfer around 
q 2 = 1 fm -2 . The inherent uncertainty when considering 
low-q data only has also been pointed out by Sick for an 
extraction of the proton charge radius from elastic elec- 
tron scattering 
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Although the value of obtained 
with our new data at low q agrees within error bars with 
the global fit it is clear from Fig. [4] that the fit of all 
available data with a suitable ansatz is more reliable. 

The deviations of the model transition densities at 
large radii (cf. Fig. [3]) lead to ( r 2 ) values exceeding 
the experimental results. Moreover, Fig. 2] demonstrates 
that also the slopes at q = (proportional to ( r A V ) are 
too steep indicating the need for an improved description 
of the tails of the model wave functions. 

To summarize, the pair decay width of the Hoyle state 
has been extracted from a new 12 C(e,e') experiment at 
low momentum transfers and independently with a novel 
model-independent global fit of the world data. The lat- 
ter, shown to be less prone to systematic errors, provides 
= (62.3 ± 2.0) fieV in accord with values deduced 



40 years ago but with much reduced uncertainty. Com- 
bined with an improved determination of the pair width 
branching ratio presently underway, the astrophysical re- 
action rate of the triple-a process is then known with a 
precision sufficient to quantitatively constrain the mod- 
eling of a variety of key problems in nuclear astrophysics. 

The result for the pair width deviates from a more 
recent Fourier Bcssel (FB) analysis. The origin of this 
discrepancy is explained by shortcomings of the FB and 
the lack of data beyond q = 3.1 fm -1 . Furthermore, the 
study of the monopole transition to the Hoyle state by 
electron scattering at low momentum transfer provides 
important constraints on models by a unique test of the 
nuclear wave function at large distances where one ex- 
pects a-cluster structures. 
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